Characterizing Temporal Distinguishability of an N-Photon State 
by Generalized Photon Bunching Effect with Multi-Photon Interference 
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The complementary principle of quantum mechanics relates qualitatively the visibility of quan- 
tum interference with path indistinguishability. Here we propose a scheme of constructive quantum 
interference involving superposition between an A''-photon state and a single-photon state to char- 
acterize quantitatively the degree of temporal distinguishability of the A^-photon state. This scheme 
is based on a generalized photon bunching effect. Such a scheme can be extended to other more 
general cases. 
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I. INTRODUCTION 

The complementary principle of quantum mechanics 
was first proposed by Bohr [l| to deal with the wave- 
particle duality of quantum particles. On the one hand, it 
successfully explained the peculiar quantum behavior of 
particles in interference. On the other hand, it only pro- 
vides a qualitative description of quantum interference 
process. The problem stems from the lack of a quantita- 
tive definition of distinguishability. Efforts were made to 
find such a definition with some success 0, d, 0, Q • 

The above mentioned discussions of the complemen- 
tary principle were mostly confined in fundamental con- 
ceptual study and in interference involving only one par- 
ticle. However, recent interests on quantum information 
involve quantum interference of multiple particles fd\ , es- 
pecially in the context of linear optical quantum com- 
puting with qubits realized by photons Q- An issue thus 
arises about distinguishability among the photons that 
may degrade the quantum interference effects, leading to 
poor performance of the quantum operations. So it is 
desirable to study photon distinguishability qualitatively 
and to find its relation with multi-photon quantum in- 
terference effect. 

The first investigation of the effect of photon distin- 
guishability on multi-photon interference was performed 
by Grice and Walmsley \§\ with an analysis on a two- 
photon polarization Hong-Ou-Mandel interferometer Q . 
A more complicated four-photon case was studied by Ou 
et al [13, [iH and later by Tsujino et al with con- 

cerns about the distinguishability between two pairs of 
photons. 

Recently, the current author ^14i] made an attempt to 
generalize to the above discussion to a system of arbi- 
trary number of photons. A degree of temporal distin- 
guishability is quantitatively defined and a destructive 
multi-photon interference method is proposed that relies 
on a quantum state projection measurement ITgI . [iTj 
to experimentally measure it. Subsequent experimental 
demonstrations [isl . llOll confirmed some of the predic- 
tions. It was shown [14l| that the visibility of interference 



is proportional the number of indistinguishable photons 
in a simple situation. But since the visibility is bounded 
by 1, more accurate measurement on the visibility is 
required to distinct various scenarios of different pho- 
ton distributions, especially when the photon number is 
large. The accuracy problem is compounded by the fact 
that destructive interference in this scheme makes the 
measured quantity small, so that it requires long record- 
ing time for good accuracy. Furthermore, the scheme of 
quantum state projection measurement is complicated in 
structure and requires phase shifters with precise values. 
Another scheme was recently discussed in Ref.flOl that 
relies on a generalized Hong-Ou-Mandel interference ef- 
fect with asymmetric beam splitter. This scheme needs 
less optical elements and thus significantly simplifies the 
optical arrangement. But since this new scheme is based 
on destructive interference, it still suffers the problem of 
low count rate at maximum interference effect and thus 
low measurement accuracy. Another disadvantage in this 
new scheme is that we need to control the precise value 
of the transmissivity of the beam splitter, which depends 
on the total photon number, in order to achieve complete 
destructive multi-photon interference. 

It was demonstrated that a photon pair bunching effect 
can be used to characterize t emp oral distinguishability 
between two pairs of photons [lOillH. In this case, con- 
structive four-photon interference for two pairs of pho- 
tons leads to five-fold increase in four-photon coincidence 
when the two pairs are indistinguishable, whereas the 
state of two separated pairs produces only three-fold in- 
crease. The enhancement factor is expected to be bigger 
for larger photon number due to Bose statistics. Large- 
ness of the measured quantity leads to a good accuracy 
in the four-photon coincidence measurement. 

In this paper, we generalize the study of the photon 
bunching effect to arbitrary photon number. Wc find 
that the enhancement factor in photon bunching is due 
to constructive interference and can be used to charac- 
terize the temporal distinguishability of photons. The 
enhancement factor is not sensitive to the experimental 
parameters and the optical arrangement is relatively sim- 
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pie. This scheme seems to overcome all the shortcomings 
of previous schemes. The paper is arranged as follows: In 
Sect. II, we will discuss stimulated emission as a photon 
bunching effect due to constructive interference and ex- 
ploit it for characterizing temporal distinguishability of 
incoming photons. Wc will also discuss its analogy with 
a beam splitter. This is a simple single mode analysis. In 
Sect. Ill, we will perform the more rigorous multi-mode 
analysis and confirm the results from the simple single 
mode analysis. In Sect. IV, we will consider other sce- 
narios of photon temporal distributions and derive the 
corresponding enhancement factor. We end the paper 
with a discussion. 



II. STIMULATED EMISSION AS A 
MULTI-PHOTON CONSTRUCTIVE 
INTERFERENCE EFFECT 

Recently, it was pointed out [2l| that stimulated emis- 
sion can be interpreted as a result of multi-photon con- 
structive interference: when N input photons are in- 
distinguishable from the photon emitted by the excited 
atom, constructive interference leads to a factor of N en- 
hancement in the atomic emission rate from spontaneous 
emission. The enhanced emission is due to stimulated 
emission. On the other hand, if the input photons are 
completely distinguishable from the photon emitted by 
the atom, no enhancement occurs and the atom makes 
only spontaneous emission. 

If the input photons are partially indistinguishable 
from the emitted photon, only the indistinguishable part 
will give rise to the stimulated emission. Therefore, such 
a scheme can be used to characterize quantitatively the 
degree of distinguishability of the input photons. To see 
how this works, we consider an excited atom modelled 
as a phase insensitive quantum amplifier with small gain 



gal, 



(1) 



where So represents all the internal modes of the amplifier 
and it is usually independent of the signal mode and 
is in vacuum. To preserve the commutation relation, we 
need |Gp — |gp = 1 and for small gain, \g\ « 1. The 
related evolution operator for the system has the form of 



U = exp{f]alal — h.c.} w 1 + (c/ajaj + h.c.) (2) 
with g ~ T]. 

With a vacuum input of |0), we have the output state 



t;;t 



^)'^l = U\0)^\0)+9\l)s®\l)o- 



(3) 



This gives the spontaneous emission probability of \g\'^. 
When the input is an A^-photon state |iV)s(g) |0)o, we have 

\^)^ol « \N),\0)o + gial\N),)®{al\0)o) 

= \N),\0)o + gVNTl\N + l),®\l)o. (4) 



The probability becomes {N + l)\g\'^- The stimulated 
emission helps to enhance the emission rate by a factor 
of AT + 1. 

In Eq.(|3]), the input N photons are all in the same 
mode as the mode of the amplifier. However, if some 
of the input photons are in different modes from the 
mode cis of the amplifier, these photons are not cou- 
pled to the amplifier and cannot stimulate the emission 
of the amplifier. Mathematically, we have the input as 
|m)s|A^ — to)s'|0)o and the output state as 



\m),\N ~m)s'\0)Q 

+g{al\m)s) (E) {\N - m) ® (aj|0)o) 
= |m),|Ar-m)£'|0)o 

4-5Vm + l|m + \) s\N - m),.|l)o. (5) 

The enhancement factor is now m -I- 1. In the special 
cases when m = 0,N, we recover Eqs.([31Hl), respectively. 
Therefore, spontaneous emission corresponds to the case 
when the input photons are completely distinguishable 
from the photon emitted from the amplifier whereas stim- 
ulated emission occurs when the input photons are indis- 
tinguishable from the photon emitted by the amplifier. 

Notice that the enhancement factor m -|- 1 is linearly 
related to the number of indistinguishable photons. Thus 
by observing the size of the enhancement, we can quan- 
titatively characterize the degree of distinguishability. 
However, the mode of the amplifier is somewhat com- 
plicated, which makes this scheme hard to implement. 
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FIG. 1: Generalized photon bunching effect for characterizing 
the distinguishability of photons. 

We can circumvent this problem with linear optics. As 
discussed before, the enhancement effect in stimulated 
emission is due to photon indistinguishability and is the 
result of constructive multi-photon interference, which 
can then be mimicked by a lossless beam splitter, as 
shown in Fig.l. Ref.[2l| showed that the result in the 
scheme of Fig.l for N + 1 indistinguishable photons is 
the same as the stimulated emission process described by 
Eq.Q with an enhancement factor of iV 4- 1 compared 
to the situation when the N photons are distinguishable 
from the single photon at the other side of the beam 
splitter. 

For the case when m ^ iV, i.e., the case when some of 
the input N photons are distinguishable, we may use a 
similar input state as in Eq.®, i.e., \N — m)a' \'m)a\i)b- 
In this state, the N — m photons are distinguishable from 
the m photons and the single photon from the other side. 
Since the state |m)a|l)6 is the same as the result of stim- 
ulated emission with input state of \m) and the state 



3 



\N—m)a' has no enhancement effect, the overall enhance- 
ment factor is simply m+l, exactly the same as the case 
of stimulated emission given in Eq.©. 

The enhancement effect with a beam splitter in Fig.l 
is a generalized photon bunching effect and it can be 
similarly used for the characterizing the degree of photon 
distinguishability. 

However, the above analysis is a single-mode analysis. 
To fully prove its validity, we need a multi-mode analysis. 

III. MULTI-MODE ANALYSIS OF THE 
MULTI-PHOTON BUNCHING EFFECT 

In this multi-mode analysis, we will only concentrate 
on the temporal/spectral mode and ignore other modes 
such as spatial and polarization modes. However, the 
generalization is straightforward. 

Let us consider the scheme in Fig.l where an iV-photon 
state and a single-photon state enter a beam splitter from 
two sep arate sides (labelled as a and h) . According to 
Ref. [ij , an arbitrary TV -I- 1-photon state of a wide spec- 
tral range can be expressed for multi-mode analysis as 

|$7V,i) = N^^^^ j dLJodLJidLj2..-dijjN^{LJo;iiJi, ■■■,ojn) 
xb\LUo)a\cui)d\LU2)...a^u;N)\0), (6) 
where the normalization factor J\f is given by 

TV = y duJodLJl...duJN^*{uJo',^^l, ■■■,i^n)^ 

x^$(cc;o;P{c^i,...,c^iv}) (7) 

p 

with P as the permutation operator on the indices of 1, 2, 
TV. The sum is over all possible permutations, d and 



b represent the input modes a and b of the beam splitter, 
respectively. 

Let us consider the situation when the single photon 
from input port b overlaps temporally with m photons 
from the N input photons at input port a and the rest of 
the [N — m) photons in side a are completely distinguish- 
able in time from the (m-|- 1) photons. From Ref.[14|, we 
learn that the N -I- 1-photon wave function satisfies the 
permutation symmetry relations: 

$(wo;a;i, ...,ujn) 

= $(P{a)o;'^i, ...,t^m},t^m+l, ••■,t^Af), (8) 

for all permutation operation P and the orthogonal rela- 
tions: 

j dijjQduJi...dwN^* {ujQ]uJi^ ...^LOisf) 

x$(Pfc,{wo;wi,...,c^A'}) = 0, (9) 

where P^j interchanges the indices k,j with k < m, j > 
m + 1. Here Eq.® is for indistinguishability among the 
m + 1 photons whereas Eq.® is for temporal distin- 
guishability between the m+l photons and the remaining 
N — m photons. Eqs.([51 ^ can be equivalently written 
in time domain as 

G{to; ti, ijv) 

= G{P{to; h, i™}, ijv) (10) 

and 

J dtodti...dt]yG*{to;ti, ...jt^) 

xG{Fkj{to-:ti,...,tN})^0, (11) 

where the notations are the same as before and 



Gito;h,...,tN) = (3^yjv+i)/2 J rf^o...dc^JV<i>(^o;wi,...,c^iv)e-'("«*''+"^*^+-+""*"). (12) 

The N+ 1-photon coincidence rate of the iV+ 1 detectors in Fig.l is proportional to a time integral of the correlation 
function of [2J] 

T(^+'Hto,h,...,t^) = {^^^,\E[''^\tr,)...E[°^^^ (13) 

where 

E[''\t) ^ VTEait) + VREb{t) with Ec{t) ^ {1/V2^) J duj c{Lu)e-"^* {c = a,b). (14) 
Let us first evaluate E[°\tQ)E[°\ti)...E[°\t]^)\^]\j^i) , which has the form of 

N 

E[°\t„)E[°\h)...Ei°\tN)\<i'NA)^T^/^R'/^Y.^ok^^^^ (15) 

k=0 
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It is straightforward to show that for the state in Eq.®, we have 



.+UlNtN) 



|0) 



where 



= AA-l/2^(to;il,...,f^)|0), 

g{to;ti, ...,tN) = J2G{to;nti, -^tN}). 



k=j 



(16) 
(17) 



The overah (A'' + l)-photon coincidence probabihty is then a time integral of the F- function in Eq. (fT3| : 

Pn+i= J dtodh...dtNT<^''+'Hto,h,...,tN). (18) 
With Eqs.dia [m Uni nil), we obtain 

PN+i = T^m-^ f dtodh...dtNY,Pok{g*{to,h,...,tN)}Poj{G{to,ti,...,tN)} = T^R(j2+Y.)- (^9) 

It is straightforward to find that the first term in Eq. (fT9|) is 

^=AA-l(iV + l) j dtodti...dtNYok{Q{h,tu ...,tN)}^ ^ M-\N + I) j dtQdti...dtN\Q{to,h,...,tNi\ , (20) 



where we switched the integral variables and t^. Using Eq. (|17p and the fact that is over all permutations in 
{ii, tAf}, Eq. ([20|l becomes 

^^M-\N + 1) I dtadti...dtNN\G*{to-ti.-.tN)^G{to;f{ti,...,tN})^M-^{N + l)NW^ (A^+1)!, (21) 

k=j p 

where for an arbitrary permutation in Q* ^ we also permute the variables of integral in the same way. Since Q is 
unchanged under any such permutation, we obtain A^! identical terms in Ea. (|2ip . 

Next, let us consider the second term in Eq. (ll9p . We evaluate one arbitrary term in the sum (fc ^ j): 



I dto...dtN^G*{tk,^{ti,...M,-,tN))^G{t,y{ti,...M.-.tN}) 



(22) 



Because of the permutation properties in Eqs. fTUl [TT|) 
and tk 7^ tj, the only way to obtain a non-zero integral 
is for tk in P'{ti, to, iAr} to be permuted to the first 
m positions by P'. Then we can use the permutation 
relation in Eq. pil)) to interchange it with tj so that for 
these P's, we have: 

G{tj,¥'{ti, ...,tQ, ...,t^}) = G{tk,V{ti, to, iw})- 

(23) 

Permutation by P' to other positions for tk cannot be 
interchanged with tj and by the orthogonal relation in 
Eg. pTjl . the integral is zero. Since this is only about tk 
in P', other — 1 time variables in P' are free to move. 
So, there will be m{N — 1)! permutation terms in the 
sum over P' that are nonzero and, as before in Ea. (pT|l . 
they will all have the same time integral of {N + l)J\f. 



Therefore, the second term in Eg. p^ is equal to 

'^=Af-^"^m{N -l)\Af = miN +1)1. (24) 
k^j k^j 

Combining Eqs. ([211 [24)1 . we have 

PM+i^T^R{l + m){N + (25) 

The case of m = corresponds to the situation when the 
single photon of port h is completely distinguishable from 
all the N photons from port a, which gives no interfer- 
ence and sets a baseline for reference. Thus the enhance- 
ment factor for m 7^ is m + 1, in agreement with the 
single mode analysis. Note that the enhancement factor 
does not depend on the transmissivity T but the detec- 
tion probability Pjv+i does. From Eq. p^l) . we find the 
maximum detection probability at T = N/ [N +1). 



For the intermediate case when there are some partial 
indistinguishabiHty among the photons, the wave func- 
tion G does not satisfy Eqs. lfTOl [TT|) . We wiU not obtain 
a simple close form as Eq. ips]) . However, for a special 
case when all the N photons from port a are indistin- 
guishable from each other and they are only partially 
indistinguishable from the single photon input from port 
h, we have the permutation symmetry relation: 

G{to-M,...,tN) = G(<o;P{ii,...,ijv}) (26) 

but no orthogonal relation similar to Eq. pTjl . Then 
Eq. (fT7|) becomes 

g{W,ti, ...,tN) = N\G{t^-M, ...,tN). (27) 
With some manipulations, we find Eg. p^ becomes 

Pn+1 = T^R{N + \)\{1 + NVn) 

= P^'+i(l + 7VV^) (28) 

with 

/ dtQdti...dtj\jG*{tQ; ti, tN)G{ti;to, iw) 



(a) 



J dtadti...dtN\G{to;ti, ...,tN)\'^ 
J dujo...dujN^*{uJo-uji, ...,wjv)^'(wi;wo, ww) 



/ dLOo...du;N\^{tOo]LOl,--,UJNW 



Thus the enhancement factor is 

Pn+1 



TDcl 



1 + NVn. 



(29) 
(30) 



Note that when Vat = 1, we have the maximum enhance- 
ment factor of iV -|- 1 , indicating complete indistinguisha- 
bility whereas when Vn = 0, there is no enhancement 
effect, due to complete distinguishability. So the quan- 
tity Vat gives the degree of distinguishability between the 
single photon in port b and the N photons in port a. 

The more general case is when there are m indistin- 
guishable photons among the N photons in port a and 
the m photons have partial indistinguishability from the 
single photon in port b. To simplify the case, we further 
assume that the other N — m photons from port a are 
completely distinguishable from the above m + 1 pho- 
tons. For this case, we can show similar to Eg. ((30)) that 
the enhancement factor is 

Pn+1 



pel 
^N+1 



^1 + mV„ 



(31) 



where V„i is defined in Eq. (|29p but with the wave func- 
tion G satisfying the permutation symmetry relation in 
Eq.lUni) in stead of Eq.®. 

Now the experimental procedure to measure the distin- 
guishability of the iV-photon state is depicted in Fig. 2, 
where we scan the relative delay of the single photon 
in port b with respect to the iV-photon state in port a. 
Whenever the single photon scans through m indistin- 
guishable photons, the TV -I- 1 coincidence count shows 
a bump of size m relative to the baseline. In this way, 
we can characterize the temporal distinguishability of the 
A'^-photon state. The width of the bump is determined 
by the function Vm- 



(b) 



A 



m+l 



FIG. 2: (a) A temporal distribution with well-separated 
groups of photons and (b) the corresponding normalized 
Pn+1 as the position of the single photon is scanned. 



IV. MORE GENERAL CASE OF \Na,Mb) 

For the more general input state of \Na,Mii), we can 
use the analysis similar to Eqs.ljUISl). First, if the N + M 
photons are indistinguishable, then from the quantum 
theory of a lossless beam splitter [23|], we may find the 
probability of finding all N + M photons in one output 
side of the beam splitter as 



Pn+m = T'^R'^'iN + M)\/N\M\. 



(32) 



But when the incoming N photons are distinguishable 
from the M photons, the N + M photons act like clas- 
sical particles and follow the probability law. The corre- 
sponding classical probability is then 



pel _ rpN jyM 
^N+M — ^ ^ ■ 



(33) 



So the enhancement factor due to quantum interference 



Pn+m _{N + M)\ 



pel 

^N+M 



N\M\ 



(34) 



A special case is for = Af = 2, which gives the ratio of 
6. This is the photon pair bunching effect experimentally 
demonstrated by Ou et al p^ . 

The photon bunching enhancement factor in Eq. (j34|) 
is for all the photons involved to be indistinguishable. 
When some of the photons are distinguishable, the en- 
hancement factor will decrease. The most general sce- 
nario is that some of the N photons at input port a are 
indistinguishable from some of the M photons at side b. 
Let's break the N photons and the M photons into fc -I- 1 
groups, respectively, namely, N = ni + ... + Uk + n^+i 
and RI = mi + ... + + mk+i. In these groups, the 
Hi photons are indistinguishable from rrii photons with 
i — 1,2,..., k and {ni,mi} group of photons are distin- 
guishable from {nj,mj} group of photons with i ^ j. 
Furthermore, the n-k+i photons are distinguishable from 
nik+i photons. Such a situation is depicted in Fig. 3. 

In analogy to the case of stimulated emission described 
in Eq. ^ , we may write the input state to the beam split- 
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FIG. 3: Temporal distributions for photons from input sides 
a and b, respectively. 



TABLE I: Enhancement factor for 2 a-photons and 2 b- 
photons input 



2a26 


2al6 + lb 


lab+lab 


lab+a + b 


6 


3 


4 


2 



ter as 

k 
i=l 

(35) 

where we use (8) and superscript "(*)" to separate and 
label the states of distinguishable photons. 

Since |ni)a|mi)b is the same state as \Na,Mi,) that 
gives rise to the enhancement factor in Ea. lp4)) . it will 
contribute a factor of (n^ + mi)l/ni\mi\ to the overall 
enhancement factor, which is then 

pel LL ■ 

Note that since \nk+i)a and |TOfc+i)f, are distinguishable 
states, they have no contribution to the enhancement 
factor. 

In Tables I-III, we list the enhancement factors for var- 
ious scenarios of the input states |2a,2b), |3a,2b), |3a,3;,), 
respectively. These tables are in contrast with the visi- 
bility tables in Ref.[ll]. 

V. SUMMARY AND DISCUSSION 

In this paper, we discussed a generalized photon 
bunching effect which may involve arbitrary number of 
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